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r^ , Abstract. We study the dynamics of the two-point statistics of the Kraichnan 

r^ I ensemble which describes the transport of a passive pollutant by a stochastic turbulent 

■ flow characterized by scale invariant structure functions. The fundamental equation 

"j^ , of this problem consists in the Fokker-Planck equation for the two-point correlation 

'tr. ' function of the density of particles performing spatially correlated Brownian motions 

with scale invariant correlations. This problem is equivalent to the stochastic motion of 
an effective particle driven by a generic multiplicative noise. In this paper we propose 
an alternative and more intuitive approach to the problem than the original one [1] 
i-rt ■ leading to the same conclusions. The general features of this new approach make 

^ ' possible to fit it to other more complex contexts. 
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1. Introduction 

Correlated random walks have always received interest from different scientific 
disciplines thanks to their ability to generate complex patterns and collective behaviors 
[llISlllllHinilTllH]. A case of particular relevance amounts to considering the case 
of a particle distribution in which single particles perform an ordinary Brownian 
motion, but the trajectories of different particles are spatially correlated [9l [10]. A 
rich phenomenology occurs in this case ranging from simple diffusion, to clustering and 
coalescence of particles [121 [H]- The so called Kraichnan ensemble [13], HI], introduced 
to describe the transport of passive scalars by turbulent flows, belongs to this class of 
models. In this approach, the turbulent flow is modeled as a stochastic velocity field 
with no time correlations (i.e., white in time) and scale invariant spatial correlations. 
Passive particles (e.g. pollutant) are considered to have vanishing mass (no inertia), and 
are simply transported by the flow without affecting it. As we discuss below, from the 
point of view of non-equilibrium statistical physics, this ensemble of models is equivalent 
to the general problem of a (i— dimensional Langevin equation driven by multiplicative 
noise only. The whole ensemble has been solved in [H [15] using the mathematical theory 
of boundary conditions of elliptic differential operators developed in rigorous quantum 
mechanics |16j . 

In this paper, besides showing the strict relationship between Kraichnan turbulence 
and the physics of multiplicative noise, we present an alternative approach to these 
models leading to the same quantitative classification of the phases occurring in the 
Kraichnan ensemble. Our approach deals directly with the equation for the two-point 
correlation function of the density of pollutant instead of the correlation function of the 
passive scalar. The two equations are known to be equivalent being one the adjoint of 
the other. The reason to propose this alternative solution is twofold: (i) it represents a 
more intuitive technique based on a simple regularization of the equation at small scales. 
This regularization amounts to studying the integrated mass of pollutant surrounding 
a generic pollutant particle, (ii) The new method is so general that can be potentially 
applied to more complex cases which, up to now, have been solved only under particular 
conditions [TTJ [18] . 

The key point of the present method is an "exponent hunter" technique allowing 
an exact classification of particle-particle correlations at vanishing interparticle distance 
where the dynamics presents a singularity. 

The equation we study describes three different regimes for the pollutant 
concentration which depend on: space dimensionality, anisotropy and scaling properties 
of the velocity field correlations. These regimes correspond to three different boundary 
like situations for the singularity of the equation at small interparticle separation: 

(i) Phase of strong flow compressibility. There is only one possible solution of the 
problem, as different particles collide with vanishing relative velocity. Thus particles 
coalesce in finite time and no stationary state is approached by the system; 

(ii) Phase of weak compressibility. Particles never collide and the system converges 
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towards a stationary state characterized by power law or "fractal-like" density 
correlations; 

(iii) Phase of intermediate compressibility (called also sticky phase [I5]). Particles collide 
in finite time with non-zero relative velocity. Therefore both previous cases are 
possible depending on the boundary condition, absorbing or reflecting, assigned 
"by hand" at vanishing interparticle separation. 

Note that for the first two cases there is no chance to impose arbitrary boundary 
condition at vanishing interparticle separation because the solution of the fundamental 
equation is unique. 

The paper is structured as follows: in sec. [2] we briefly review the basic properties 
of the Kraichnan ensemble. In sec. |3] we discuss the method of solution of the Fokker- 
Planck like equation describing the evolution of the two point correlation function of 
the pollutant concentration. Finally sec. H] is devoted to discussion and conclusions. 

2. The Kraichnan ensemble and basic equations 

We start the discussion by introducing briefly, in the context of stochastic processes, 
the definition of the Kraichnan ensemble. The microscopic concentration (or particle 
density) of the A^ passive pollutant particles in the d— dimensional volume V is 

N 



Pi^^t) =^^[x-Xi{t)] 



i=l 



in which Xi(t) represents the position at time t of the i*^ pollutant particle transported by 
the turbulent flow. The particle distribution is supposed to be spatially homogeneous so 
that {p{x, t)) = Po = limy^oo N/V > 0. Each particle is considered to have a negligible 
mass, and consequently it will follow the trajectory defined by the Lagrangian flow: 

^ = v[xit),t], (1) 

where v{x,t) is the Gaussian stochastic velocity field representing the "synthetic" 
turbulent flow advecting the pollutant. The spatio-temporal correlation properties of 
this field are taken to be 



vjx^t) = 

v^{x,t)v,ix',t') = 5it-t')c^,{x-x') ' ^^ 

with /i, z/ = l,...,d. In other words, the velocity field is considered white in time and 
colored in space. 

We limit our study to the description of the two-particle dynamics and correlations. 
The higher order statistics has been characterized in Ref . [l\ . The relative vector distance 
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r(t) = Xi{r) — Xj{t) between any pair of particles satisfies the closed Langevin equation 
in Ito representatiorljl: 

The quantity w{t) = v[xi{t),t] — v[xj(t),t] is a Gaussian noise whose correlation 
properties are 



^(^) = Q (4) 

w^{t)w,{t') = 2S{t-t')d^4rit)] ^ ' 

where dfj,i,{r) = [c^j,/(0) — c^u{r)] is the so-called structure tensor oi the velocity field. 
Equations ([3]) and (jl]) indicate that the quantity r{t) satisfies the most general Langevin 
equation for multiplicative noise [191 ED] • 

The Kraichnan ensemble is defined by the scale invariant choice 

df,u{r) = ar^6f,u + hr^~'^r^r^ (5) 

with a and h constant such that Cfj_y{r) is a non- negative definite rank 2 tensor. This 
means that its Fourier transform (power spectrum tensor) 



c^,u{k) = / d'^rCf,^{r)e 



has non-negative eigenvalues for all k. Such a condition implies also that < .^ < 2 |10j . 
The one-dimensional case is recovered when setting d = 1 and 6 = 0. In the paper [1] 
and all related literature, the tensor d^i,{r) is written in terms of two other constants, 
A and B, such that 

a = A+{d + i~l)B 

h = i{A-B) ■ ^' 

The positive definiteness of the covariance tensor c^y{r) is guaranteed by taking 
y4, i? > 0, which is equivalent to the condition —{^a)/{d + ^ — 1) < b < ^a. These two 
constants are introduced as A = corresponds to the incompressible case where V-v = 
and i? = to the purely potential one with v = V0. The quantities S"^ = A+ {d — 1)B 
and C^ = A are proportional to (HVvp) and ((V ■ v)"^) respectively, and they satisfy 
the inequalities S^ > C^ > 0. The degree of compressibility is defined as the ratio 
V = C^/S"^ G [0, 1]. The bounds V = and V = I define a completely incompressible 
and a completely compressible flow respectively. In one dimension, S^ = C^ > and 
therefore V = 1. Clearly, in real flows Eq. ([5]) can hold only in a limited range of scales, 
as in general c^i,{r) must vanish for r — > oo, and accordingly d^y{r) must converge to 
the constant c^iy(O). 

As a consequence of Eq. ([5]), the point r = is a stationary point of Eq. ([3]). As 
shown below, it may determine, for particular values of the parameters "P, d and ^, an 

X A simple way to understand why Eq. [T] has to be interpreted in Ito rather than Stratonovich 
representation is that the corresponding Fokker-Planck equation must not contain the drift term, as 
the stochastic advecting flow is considered isotropic. 
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absorbing phase for the dynamics. We will discuss this problem in the Fokker-Planck 
formalism. 

The Langevin equation ([3]) is equivalent [21] to the d— dimensional Fokker-Planck 
equation (FPE) 

dtP{r, t|ro, to) = Y. ^l- [df.uir)Pir, t\ro, to)] , (7) 

describing the evolution of the probability density function (PDF) P(r,t|ro,to) of the 
interparticle distance r at time t given the initial value Vq at time tg. In Eq. (^^ dt = -^ 
and 9^^ = g^^g^ with r^ being the yU*'^ component of vector r. In turbulence literature 
(e.g., see [HI [1]) instead of Eq. i^^, one usually studies the adjoint equation which is 
called, in the stochastic processes context, a "backward" FPE, while Eq. ([7]) is called the 
"forward" FPE. They are practically equivalent, even though the differential operators 
in the two equations are defined on different functional spaces [22]. The two equations 
coincide for incompressible flows. 

Equation ([7]) is a d-dimensional FPE without drift, typical of generalized 
multiplicative noise. If the tensor d^v{f) never vanishes in the interior of the domain 
of definition of Eq. ([7]), a unique solution exists when the initial [e.g. P(r, to|ro, to) = 
5{r — Tq)] and the boundary conditions are given [211 [23]. However, if d^y{r) vanishes at 
some point, the solution could be no more unique and an additional boundary condition 
might be necessary at such a singularity. 

We develop a novel method to provide, not only the qualitative classification of 



the singularity r = for the Kraichnan ensemble (see Appendix A ), but also the exact 
solution of Eq. (I7j) in the form of a particular series expansion in the neighborhood of 
the singularity. This represents an alternative and more intuitive approach than the 
original one [1]. This method can be potentially generalized and applied to different 
kind of FPE with singular points in the domain of definition. 

It can happen that the singularity r = behaves as an additional regular boundary 
Then an absorbing, reflecting or mixed boundary condition needs to be assigned. In 
all the other cases the singularity can be seen as a boundary at which the condition is 
unique and automatically fixed by the form of Eq. ([7]) itself. In particular when the 



singularity corresponds to an adhesive boundary (see Appendix A ), the diffusion of the 
effective particle described by Eq. ([7]) is such that the boundary is reached in a finite 
time and it behaves automatically as an absorbing boundary [23]. Naively one can say 
that the effective particle, whose position is r{t), arrives at this boundary in a finite time 
but with vanishing velocity obliging it to stay at that point indefinitely. Instead when 
the singularity behaves as a natural attractive boundary the same effective particle 

§ We use here the definitions given by Van Kampen in Chap. XII of its book [23]. This classification 
differs from that given originaUy by Feller |22] and re- proposed by Gardiner [21]. The two are related 
directly by the following simple positions: "adhesive boundary" becomes "exit boundary" in the latter, 
"natural repulsive" becomes "entrance" or "natural" depending if the average time to be expelled at 
finite distance from the singularity is respectively finite or infinite, and finally "natural attractive" also 
becomes simply "natural" . 
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approaches it, but only in an infinite time. Therefore there is no need of external 
boundary conditions at the singularity. Finally, for a singularity which is a natural 
repulsive boundary for Eq. ([71), the effective particle, when placed arbitrarily close to 
this point, is expelled far from it. We will see that this last situation, in turn, splits 
into two sub-cases in which the time to be expelled at a finite distance is respectively 
finite or infinite. The details about the Van Kampen classification will be discussed in 



Appendix A 



3. Method and solution 

In this paper, instead of considering the PDF P(r, tl^o, to); we study the equivalent 
problem of the evolution of the two-point correlation function of the pollutant particle 
density at time t 

r(r,t) = MM^(^±M)_5(r), (8) 

Po 

where we have subtracted the Dirac delta contribution in r = 0, to exclude the self- 
correlation of each particle with itself [23]. In this way T{r,t) gives exactly the average 
conditional density of other particles seen by a generic particle at separation r from 
itself. The symbol (...) indicates the ensemble average over the realizations of the process 
(the hypothesis of statistical spatial homogeneity is done implicitly). It is noteworthy 
that this average in Eq. ([8]) implies both the average (...) over the realizations of the 
velocity field and the ensemble average over the initial particle configurations. From the 
knowledge of P(r, t|ro, to), we can write: 

r(r, t)= j (fro r(ro, to)P{r, t|ro, to) • 

Consequently, r(r, t) satisfies the same forward FPE as P(r, t|ro, to); which can be cast 
to the tensorial form: 

dtV{r,t) = V ■ [V ■[d{r)V{r,t)]] , (9) 

where d{r) is the flow structure tensor which, in intrinsic form, reads 

d{r) =ar^i + hr^''^{r®r) (10) 

with / being the identity tensor and r ® r the tensor product of r for itself. Let us 
now integrate explicitly Eq. ([9]) over the sphere S{r) of radius r around the origin. The 
quantity 

N{r,t)= f d\'T{r',t)= f dr' r"^~^ f dnT{r',t) (11) 

Js{r) Jo Jut 

is the number of other particles seen in average by a particle in the sphere of radius r 
centered on it, and Q^ is the total solid angle in d dimensions. Note that from Eq. flTTl) 
we have 

dnT{r,t) = ^drN{r,t). (12) 
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Using Eq. (TTT1) and the divergence theorem of integral calculus, we can write Eq. ^ as 

dtN{T,t) = f dS ■ {V ■ [d{r)T{r,t)]] , (13) 

JdS{r) *- -" 

where dS{r) indicates the surface of the sphere S{r), and dS = r'^^^dflir is the vector 
surface element of the sphere of radius r along the radial unitary vector i^. 

In order to transform the integral in Eq. ( fT3l) with the condition (ITOl) into a more 
treatable expression, we apply the tensor calculus in d = 3, and then we see how to 
extend the description to arbitrary dimensions d. The three-dimensional vector gradient 
operator in spherical coordinates in the local orthogonal spherical reference frame, is 

V = irdr + -de H ^d^ . 

r r sm & 

In the same frame, the tensor product r ® t acquires the simple form 

r 0r = r^(ir ® ir) , 
while the identity operator is obviously 

I = ir®lr + le^h + ^®^■ 
Note that i.r(e,0) ® ir(e,(/i) is the local projector operator along the direction ir(9,<f>)- 

Using the expression flTUl) and the rules of tensor calculus in the spherical 
coordinates frame [jjl, one can write 

V ■ [dir)T{r, t)] = [{a + 6)9^(r«r) + 2br^-^T] i^ 
[OeTjie 



+ ar^-' {deT)ig + ^l^ , (14) 

and therefore simply 

i, ■ {V-[d{r)T{r,t)]} = {a + b)dr[r^T{r,t)] + 2br^-^T{r,t) . (15) 

Note that the radial direction is the only important spherical component of Eq. flT4|) . 
This is exactly the reason why we have used the local spherical orthogonal frame. 
In d dimensions the gradient operator in local hyper-spherical coordinates is 
d-i 
V = Irdr + ^ ivj V,/,, , (16) 

i=i 
where V^ . is the component of the gradient along the orthogonal angular direction i^ . on 
the d— dimensional sphere of radius r. Since the curvature of the sphere is constant and 
the hyper-spherical frame is orthogonal, it is simple to show that for any j = 1, ..., rf — 1 

II In particular we use that any derivatives of the identity operator / trivially vanishes, that drir — 0, 
dgir = ig and d^ir = sinSi^, and that ir,e,4> <8) ir,e,4> are projection operators. 
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\7 T - ^'^^^ 

V,/,Jr — 



and therefore 

This observation permits to generahze the resuh flTSj) to any dimension d in the 
following way: 

i, ■ {V ■ Kr)r(r, t)]} = {a + b)dr[r^T{r, t)] + {d - l)br^-'T{r, t) . (17) 

In this relation, the only dependence on the angular variables lies on r(r,t). Plugging 
Eqs. ( TTTIl and (TT2|) into Eq. (fT3|) . we finally arrive at the following closed equation for 
N{r,ty. 

dtN{r,t) = {a + by-^dr [r^-'^+^drN{r,t)] 

+ {d-l)br^-^drN{r,t). (18) 

The singularity r = of Eq. IQ now constitutes the boundary r = of Eq. (TTSj) (the 
other one being r —>■ cxo). Equation (ITSl) can be rewritten in a form which is useful for 
the singularity classification following Feller [22] or Van Kampen [23j ( Appendix A] ), and 



to face explicitly the limit case ^ = 2 ( Appendix B ). Let us introduce n{r, t) = drN{r, t) 



such that n{r,t)dr is the average number of particles seen by the particle at the origin 
in the spherical shell of radius r and infinitesimal thickness dr around it. Equation (TT8l) 
in terms of this function becomes a one-dimensional Fokker- Planck equation 

dtn{r,t) = {a + b)d'^[r^n{r,t)] - {d - l)adr[r^-^n{r,t)] , (19) 

which depends explicitly on the dimension d only through the coefficient {d — 1) in the 
right hand side. 

3.1. The exponent hunter method for the "rough" case ^ < 2 

We now introduce a new technique for the rough fiow case .^ < 2 through which we 
not only recover the above Van Kampen's classification, but also find the exact small r 
behavior of all the possible solutions to Eq. (TT8|) . and therefore to Eq. (Q. The smooth 



case ^ = 2 is presented in Appendix B 



Since N[r,t) is by definition smooth and finite at small r, this method basically 
amounts to looking for the solution of Eq. flTSj) in the form of the appropriate power 
expansion for N{r, t). Let us assume ^ < 2 and that N{r, t) at any time can be expanded 
at small r in some power series: 

oo 

N{r, t) = Y^ ciityf^' with (3i < (3j if i < j , (20) 

1=0 

^ Due to the constant curvature of the sphere and the orthogonahty of the coordinate frame, in order 
to find V^ ir , it is sufficient to consider the equatorial circumference in the direction i^ . and therefore 
the variation of i^ along this. Clearly, as all equatorial circumferences passing through a point of the 
spherical surface are a rotation one of each other, the modulus of V^ ir is the same for any direction 
i^. and depends only on r. It is finally simple to find that such a dependence is 1/r. 
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where all f3i are independent oft. Then we construct a recursive solution, by substituting 
this expression into Eq. ([Ti 



1=0 1=0 

+ id-i)b]ci{ty+'''-\ (21) 

The only possibility to have a solution of Eq. ([21]), being ^ < 2, is that either (i) /3o = 
or (ii) [{a + b){^ + (3o — d) + {d— l)b] = 00 In order to accept or not these conditions, we 
have to take into account the statistical and physical meaning of N{r,t). This requires 
N{r, t) has to be finite and non-negative for all finite r and t. Moreover drN{r, t) has to 
be the same non-negative, as the number of particles seen cannot decrease by increasing 
the radius of the sphere. Let us analyze separately the two cases: 

(i) Assume /?o = 0. The series (120|) satisfying Eq. (12T1) needs to have Vn > coefficients 
such that: 

r Pn = n{2 - 



< 



c„(t) = 7nC„+i(t) ^22) 



7n= (n + l)(2-0{(a + 6)[^-rf+(n + l)(2-0] 
+ (rf-l)6} 

This solution corresponds to an expansion into integer powers of r^~'^. Moreover, we 
see immediately that, as /3j < Pj for z < j, /?o = implies N{r, t) to be always finite 
at finite r and t. Note also that, as (a + 6) > and ^ < 2, the inequality 7, < 7^ 
holds for i < j too. Therefore the necessary and sufficient condition preventing 
both N{r, t) and drN{r, t) from becoming negative at sufficiently small r is 70 > 0, 
implying that 7i > for all integer i. Hence, the inequality 

70 = (2 - 0[(a + b){-d + 2) + {d- l)b] > 

provides the only condition of acceptability of this solution. It can be rewritten in 
terms of the positive ratio a /{a + b) as 

a I , ^ 

< :^^ > (23) 



a + b ~ d — 1 

i.e., V > {d ~ 2 + ^)/(2^) for the compressibility. The most important feature of 
the solution with j3q = is that, when V > {d — 2 + C,)/{2C,) strictly, the function 
N{r,t) converges for r — i> to a positive value Co(t) growing with time. This 
means that particles coalesce more and more in time forming massive point clusters. 
Indeed the condition co(t) > means that the average conditional particle density 
r{r,t) develops a contribution proportional to the Dirac delta function S{r) whose 

^ Otherwise the lowest order power at the right-hand side of the equation cannot be matched by any 
term of the left-hand side 
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increasing weight is exactly Co(t). From the point of view of diffusion theory this 
corresponds either to an adhesive or to a regular absorbing boundary at r = for 
Eq. ( IT9l) . The second important feature is given by the first correction ci{t)r^^, 
with /3i = (2 — ^), to the leading constant term of N{r,t). This implies that at 
sufficiently small r > the functioio T(r,t) is proportional to r^~^~'^ with the 
amplitude proportional to Ci(t). 

Instead for P = (rf — 2 + ^)/(2^) we have 7o = identically. As a consequence, 
if co(0) = it will be co(t) = and A^(0,t) = at all t. That is T{r,t) develops 
no Dirac delta function contribution and no coalescence occurs. In other words 
the probability to find any two particles at vanishing distance is zero at all time. 
However, at small r the function r(r,t) is proportional to r^~^~'^ (i.e. diverging 
and integrable) which denotes a simple clustering of the particle distribution. This 
behavior is typical of the singularity r = acting as either a repulsive entrance 
boundary, or a refiecting regular boundary. 

(ii) Assume now that (a + b){^ + Po — d) + {d— 1)6 = 0. The related solution of Eq. (pT!) 
is given by the equalities 



f P J C id~i)b 


= 1- 


-e + - 


d-l)a 


Pu - a s a+b - 


a+b 


/3„ = /3o + n(2-0 









(24) 



Cn{t) — A„C„+i(t) 



An = [/3o + (n + 1)(2 - 0] [{n + 1)(2 - 0(a + b)] 

Again the solution can be found as an expansion in integer powers of r^~^, but there 
is also a "singular" multiplicative contribution r^°. In order to have finite N{r,t) 
at finite r the following restriction has to be imposed: 

a + b 
which can be recast to 

> -r^ ' (25) 



a + b ~ d ~ 1 

i.e., V < (i/^^ in terms of compressibility. This is the only necessary and sufficient 
condition for the acceptability of this solution. In fact, since ^ < 2, if /5o > we 
have from the last of Eq. ( 12^ A„ > for all n, and this guarantees both N{r, t) 
and drN{r, t) to be finite and non-negative. 

When a/(a + 6) > (1-0/(^^-1) (i-e. V <d/(^) strictly, one has A^(0, t) = for all 
t and no coalescence occurs. In other words different particles have zero probability 
at any time to be found at the same spatial point. In terms of diffusion theory this 
is possible only if the point r = behaves either as a repulsive entrance boundary 
or as a regular refiective one. Moreover at sufficiently small scale N[r,t) ~ CQ{t)r^° 

Or more precisely to its angular average (1/57^) / dfl T{r, t). 
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with Po = [d — C, — b{d — I) /{a + b)] > 0. This imphes that T{r, t) is proportional 
to Co(t)r^"~'^. Note that the case a/ {a + 6) = l/{d — 1) above is in the present class 
and in fact the exponents coincide for this value of a /{a + b). The quantity D = j3q 
plays the role of a local (density) fractal dimension of the particle distribution. 
Finally, for a/{a + b) = {l — ^)/{d — l) (i.e. V = d/^'^) we have again /5o = and this 
solution belongs to the class of solutions with the adhesive or absorbing behavior 
at r = seen above for a/{a + b) > l/{d — 1) and developing particle coalescence. 

At this point we can draw the following conclusions for the rough stochastic flows 
with ^ < 2. The physics described by Eq. ((71) can be summarized in three different 
regimes (phases) depending on the values of the compressibility parameter V: 
(i) For P < (d — 2 + ,^)/(2^) (weak compressibility) the solution is unique and coincides 
with the "repulsive" solution characterized by /3o = [d — i ~ b{d — l)/(a + b)] > 
showing no coalescence but simple clustering. This conflrms that for these values of 
the compressibility V, the singularity at r = acts as a repulsive entrance boundary: 
particles never collide. 

(ii) For (d — 2 + ,^)/(2^) <V< d/^"^ (intermediate compressibility) both attractive and 
repulsive solutions are possible, and clearly also all their compositions. This confirms 
that the singularity r = works as a regular boundary because particles hit one 
each other in finite time with non-zero relative velocity. Accordingly, it is necessary 
to select an absorbing, refiecting or mixed boundary condition to determine a single 
solution. In realistic situation, the choiche of appropriate boundary conditions is, of 
course, suggested by physical considerations on the specificity of the interactions among 
pollutant particles in the fiuid environment. Moreover, in real fiows, an important role 
is played by the existence, at small separation, of two characteristic scales Cy, £« related 
to purely viscous and diffusive motion respectively. Below them, the Kraichnan picture 
does not hold anymore. The attractive solution has to be selected when {Iv^Id} — * 
in such a way the Schmidt number Sc = zz/k diverges fast enough [151 ISSl 126] . 
(iii) Finally, for V > d/^"^ (strong compressibility) the solution is again unique showing 
particle coalescence signed by the development of an increasing delta function at r = 
in r(r,t). It confirms that for such strong compressibility, r = works as an adhesive 
(or exit) boundary for which particles collide in finite time but with vanishing relative 
velocity. 

It is important to stress that, due to the structure of the equation system for the 
coefficients c„(t), the only possible stationary solution is given by Cq > and c„ = 
for 77, > 1. Clearly this is a real and physically meaningful stationary solution only 
for the natural repulsive case at low compressibility V < (d — 2 + ^)/(2^) or for the 
regular boundary case at intermediate compressibility {d — 2 + ^)/(2S,) <V < (i^/2^ with 
refiecting boundary condition. For the other cases this stationary solution is unphysical 
amounting to placing initially all the system particles at the same spatial point. It is 
in fact known from the theory of stochastic processes [23| that when either an adhesive 
boundary or a regular absorbing boundary are present no stationary solution is possible. 
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4. Discussion and conclusions 

In this paper we have presented a different approach to analyze and classify the solutions 
of the Kraichnan ensemble. Differently from the previous methods [H [15], our approach 
focuses on the equation for the two-point correlation function of the pollutant density 
instead of the correlation function of the passive scalar which is the customary quantity 
studied in turbulence. The two equations are known to be equivalent being one the 
adjoint of the other. The reason to suggest this alternative method is two-fold. First, it 
represents a more intuitive approach to the problem based on a natural regularization 
of the basic diffusion equation at small scales. Second, the new method is so general 
that can be potentially applied to more complex flows. 

The key points of the present approach are: (i) transforming the fundamental 
diffusion equation for two-point correlations r(r,t) into an equation for the integrated 
mass of pollutant A^(r, t) surrounding a generic pollutant particle for which the 
singularity is smoother; (ii) developing an "exponent hunter" technique, consisting 
in finding the appropriate power series expansion allowing an exact and quantitative 
classification of the particle-particle correlations at finite inter-particle distance. 

In this way all the possible behaviors of the small separation singularity is obtained 
directly from the explicit solution of the equation for N{r,t), and the classification 
becomes straight forward. 

Finally, it is noteworthy to observe that the crossover from intermediate to 
strong compressibility corresponds in multiplicative noise field theories to the the non- 
equilibrium second order transition from an active to an absorbing phase. In particular 
the adhesive behavior of the singularity for strong compressibility Kraichnan models 
stands for the absorbing phase in field theories [19]. 

Appendix A. The Van Kampen's classification 

In this appendix we apply to Kraichnan ensemble the boundary (or singularity) 
classification introduced by Van Kampen in [23] for the general one-dimensional Fokker- 
Planck equation 

dtf{r,t) = ^d'AD{r)f{r,t)] - dr[K{r)f{r,t)] . (A.l) 

describing the evolution of the PDF /(r, t) of the position r of a particle at time t. Indeed 
Eq. flTOl) is exactly of this type with D{r) = 2{a + b)r^ and K{r) = {d — l)ar^~^ and 
with the singularity to be classified r = 0. Van Kampen's classification for a singularity 
r = is based on the analysis of the behavior for e — >■ of the integrals: 

L. = rdre^ir,gdr--^ (A.2) 

Ls = rdr^ 
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where 

and with x, tq > 0. In general one can show that 

(i) If for e — i> we have Li -^ +oo, the singularity r = behaves as a natural repulsive 
boundary as the particle starting from tq > has zero probability to reach r = 
|23j . This means that starting the motion close to r = 0, the particle run away from 
the singularity never touching it. Consequently, Eq. flA.ll) need not an additional 
boundary condition at r = and the solution, once the initial and the other possible 
boundary conditions are given, is unique. 

In general one distinguishes two sub-cases depending whether the conditional mean 
escape time from the singularity to a finite distance is finite or infinite. In the 
first sub-case one can show that the solution converges in time to a stationary 
state. The conditional mean escape time can be evaluated as follows. Let us put at 
r = e > (which is a regular point) a reflective boundary and start the dynamics 
from a generic tq G (e, ri) with ri > e. It is possible to show [23j that the mean 
escape time through ri conditioned to starting the dynamics from tq and having a 
reflecting boundary at e is 

r(ri,e|ro)= ^;°^^^^^(^,^ . (A.3) 

We have an entrance boundary at r = if r remains finite for r ^ e ^ 0. Instead, 
if r diverges, r = is a proper natural repulsive boundary. We show below how to 
distinguish these two possibilities in our Kraichnan case. 

(ii) If for e — !► we have Li < +cxd and L2 — > +00, the point r = behaves as a natural 
attractive boundary. The particle starting at finite r approaches r = but in an 
infinite mean time. As above the solution is unique and no boundary is needed to 
be fixed. Differently from above, no stationary state is reached. 

(iii) If for e — ^ we have Li, L2 < +cxd and L3 —>■ +00, the point r = behaves as an 
adhesive boundary. The particle starting at finite r approaches r = in a finite 
time but reaches it with vanishing velocity. Therefore once the particle has reached 
r = it stays there forever. In other words r = works "naturally" as an absorbing 
boundary. Again the solution is unique and no boundary condition has to be fixed 
by hand. Since the probability to find at t > the particle at r = is finite and 
increases in time, no stationary state is reached and f{r,t) develops a Dirac delta 
function at r = with a time increasing coefficient; 

(iv) Finally, if in the same limit Li, L2, L3 < +cxo, the point r = behaves as a regular 
boundary. One can now show that the particle reaches r = in a finite time but 
with a non-zero velocity. Consequently, a solution to the equation is determined 
once a boundary condition at r = is explicitly fixed. This condition can be either 
absorbing, or refiecting, or mixed. Only if a purely reflecting condition is fixed the 
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solution runs towards a stationary state. Otherwise, as in the previous case, a time 
increasing Dirac delta contribution appears at r = 0. 

Before applying this classification to the Kraichnan ensemble, it is useful to remind 
again that in this case D{r) = 2{a + b)r^, K{r) = [d— l)ar^~^ and r is not the position 
of a single particle but the relative distance between two system particles, so that r = 
means a collision between these two particles. The above classification can be translated 
as follows: 

Case 1 is obtained for a/{a + b) > l/{d - 1), i.e., ioi P < {d - 2 + 0/(20 (weak 
compressibility). Moreover it is simple to show that for ,^ < 2 (rough fiow) the 
conditional mean escape time from r = 0, Eq. ( JA.31) . for ro — > e — > is finite and it 
behaves as an entrance boundary, while for ^ = 2 (smooth fiow) this time diverges and 
r = behaves as a proper natural repulsive boundary. In other words for ^ < 2 nearby 
particles almost surely never collide and get far away from one each other in finite time, 
while for ^ = 2 nearby particles again never collide but increase their relative distance 
logarithmically. 

Case 2 is obtained only when ^ = 2 and a/{a + b) < l/{d — 1), i.e., V > d/i (strong 
compressibility for smooth fiow). Now any pair of particles at finite relative distance 
almost surely do not collide approaching one each other very slowly. The difference 
between the two sub-cases with ^ = 2 is very subtle. This is the reason way in the original 
Feller's classification they were included in a unique class of "natural boundaries" . We 
Appendix B|how to distinguish these solutions. 



see m 



Case 3 is obtained when .^ < 2 and a/{a + b) < (^ — l)/{d — 1), i.e. V > d/^"^ (strong 
compressibility). In this case any two particles at finite initial relative distance almost 
surely collide in a finite time, but with vanishing relative velocity and consequently 
coalesce. 

Finally case 4 is obtained for ^ < 2 and (^ — l)/{d — 1) < a/{a + b) < l/{d— 1), i.e., for 
(d — 2 + ^/(20 < V < d/^^ (intermediate compressibility). Now any pairs of particle 
almost surely collide in a finite time with non-vanishing relative velocity. Therefore to 
fix the solution of the equation one as to decide if collisions are either completely elastic 
(refiecting boundary at r = 0) or completely inelastic (absorbing boundary at r = 0) or 
intermediate (mixed boundary). This is the only case in which it is necessary to fix by 
hand a boundary condition at r = to determine the solution of Eq. ([9]). 

All this qualitative analysis coincides with the one given by [Ij through the 
boundary condition theory of elliptic operators. 

Appendix B. The "smooth" case ^ = 2 

We have seen that in Van Kampen's classification the solution of the "smooth" case 
,^ = 2 to our problem is unique and corresponds to a singularity at r = behaving as a 
natural (attractive or repulsive) boundary. In rf = 1 the problem ^ = 2 has been solved 
by [Hj which found the solution to Eq. ([7]) for tq — >■ 0. The (i— dimensional case has 
been extensively studied in [27] by analyzing the time evolution of the moments of the 
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separation between two system particles. Here we give the exact solution for the average 
conditional radial density of particles whose evolution is described by Eq. flTOl) . This 
direct solution clarifies the meaning of r = as a natural boundary and the distinction 
between an attractive and a repulsive case. 

First of all we see that any power law function ri{t)r"' is a solution of the Eq. (JTOjl 
with (^ = 2 if ri{t) satisfies: 

r/(t) = r7(0)e^(°)* (B.l) 

with 7(a) = (a + l)[(a + b){a + 2) — {d — l)a]. By the physical definition of 
n{r,t) the exponent a needs to be larger than —1 to have always a finite conditional 
number of particles N{r,t) in any sphere of finite radius. Therefore 7(a) > if 
a/{a + b) < (a + 2)/{d — 1). And in particular for a/{a + b) < l/{d — 1) we have 
7(a) > for all permitted a. This gives a first insight into the difference between 
attractive and repulsive natural boundary at r = 0. In fact for a/{a + b) > l/{d — 1) 
(repulsive case) there are initial conditions which are depleted by the dynamics, while 
for a/ (a + 6) > l/{d — 1) (attractive case) all physical initial condition are amplified. 
However the fact that any spatial power law with a coefficient satisfying Eq. fIB.ip is 
always a solution of Eq. flT9l) looks unphysical because at any time the average particle 
density in the infinite volume is not conserved. This unphysical aspect is due to the fact 
that d{r) is divergent on large separations. Therefore our model has to be interpreted 
as valid up to an upper cutoff L. 

In order to better clarify the nature of the solution of Eq. flT9|) we therefore study 
the evolution of n{r, t) with initial condition 

{7r° for r < Tc 
(B.2) 
for r > Tc 

where a > — 1 and r^ is an arbitrary finite scale. 

Equation (TT9l) can be transformed into a more treatable form performing the change 
of variable r = e" and considering the function m{u,t) defined by the relation 

m(u, t) = n{r = e", t)^!^ = n{r = e", t)e" , 

du 

which conserves the measure. Directly from Eq. fITI?]) it is simple to show that 

m{u, t) satisfies the simple Fokker-Planck equation at constant coefficients typical of 

ordinary Brownian motion with diffusion coefficient 2{a + b) and constant drift velocity 

-[b+{2-d)a] 

dtm{u, t) = {a + b)dlm{u, t) + [b+ {2 — d)a]dum{u, t) . (B.3) 

The initial condition is given by changing variables in Eq. (IB. 21) : 

^g(a+l)u ^Qj. u < Uq = log Tc 

m{u,0) = { (B.4) 

ioi U > Uq = logTc 
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The solution of Eq. (1B.3P is easily found by considering the Fourier transform m(g, t) = 
J_°° dum{u,t)e~'''^'^ which, substituted in Eq. (IB. 30 . reads 

m(g, t) = m{q, 0) exp{-[(a + b)q^ - i{h + (2 - d)a)q]t} . (B.5) 

From Eq. ( IB. 41) . m(g, 0) is given by 



m(g, 0) = 7 



g(«+l-J(j)MO 



a + \ — iq 
Plugging this expression into Eq. (IB. 50 and inverting the Fourier transform one finds: 

m('U, t) = 7exp[(a + l)(-u + Cot)] / e~^ (B.6) 

JvUu) V27r 



where 



vm{u) V27r 

Co = 6 + (2 - c?)a + (a + l)(a + b) 

u-Uo + [b+{2-d)a + 2{a + l){a + b)]t 
v/2(a + b)t 
We now study the asymptotic of Eq. (IB. 61) by using the following well known 
approximat ions : 



,.2 



+00 ^^ „2 I -7s=— for Wm > 1 



/ 27ri)„ 

.g 2 ~ 



^27r 

for t>m ^ — 1 

Using this and moving back from m(u, t) to ra(r, t), we finally find for r/r^ ^ exp[— C2t + 
v/2(a + b)t] 

(log-^+Cit)2- 



;a + b)t (Tc 



exp 



n(r,t)^7r."W^^^-^(^) ,^ /!"^t ^ (B-T) 



TT V r / log f + C2t 



Tc 



with 



ci = 6 + (2 - ci)a 



C2 = 6 + (2 - c?)a + 2(a + l)(a + 6) , 

Eq. ( ]B.7p practically says that, in the region of validity of the approximation, n{r,t) 
develops a log-normal behavior whose peak drifts with velocity —Cit. 
Instead for r/r^ <^ exp[— C2t — ^/2{a + b)t] we have simply 

n{r,t) ~ 7r"exp{(a + l)cot} = n(r, 0) exp{(Q; + l)cot} , (B-8) 

i.e., the same amplifying aforementioned behavior for the scale invariant initial 
condition. 

By looking at Eq. (IB. 71) we can appreciate better the meaning of "attractive" or 
"repulsive" natural boundary behavior at the singularity r = 0. 

(i) For Ci > 0, i.e., V > d/4, the peak of the log-normal function shifts to smaller and 
smaller scales denoting the attractive behavior of the point r = 0; 
(ii) for Ci < 0, i.e, V < d/A, such peaks moves away from the singularity signing the 
repulsive behavior of such singularity. 
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